Pairing instabilities in quasi-two-dimensional Fermi gases 
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We study non-equilibrium dynamics of ultracold two-component Fermi gases in low-dimensional geometries 
after the interactions are quenclied from a weakly interacting to a strongly interacting regime. We develop 
a T-matrix formalism that takes into account the interplay between Pauli blocking and tight confinement in 
low-dimensional geometries. We employ our formalism to study the formation of molecules in quasi-two- 
dimensional Fermi gases near Feshbach resonance and show that the rate at which molecules form depends 
strongly on the transverse confinement. Furthermore, Pauli blocking gives rise to a sizable correction to the 
binding energy of molecules. 
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I. INTRODUCTION 



Theoretical prediction and experimental observation of 
magnetic-field-induced Feshbach resonances in ultracold 
atoms [1] paved the way for many exciting discoveries, in- 
cluding demonstration of fermionic superfluidiU' iQ], observa- 
tion of Efimov trimers and Fermi polarons [3h31, and creation 
of quantum degenerate gases of polar molecules Surpris- 
ingly, additional Feshbach resonances can be found in sys- 
tems with reduced dimensionality. Earlier theoretical work on 
two particle scattering in systems confined to one-dimensional 
tubes l'9'l and two-dimensional pancakes [lO] suggested a pos- 
sibility realizing "confinement-induced resonances" (CIR's), 
i.e., special scattering resonances made possible by restrict- 
ing the transverse motion of atoms. Such resonances have 
been observed in both one-dimensional (ID) ifTlll and two- 
dimensional (2D) il2ll systems. In optical lattices, Feshbach 
resonances can give rise to nontrivial manifestations of mixing 
of higher Bloch bands iflsilTl. 

Most of the earlier work has focused on the interplay of 
dimensional confinement and resonant interactions in two- 
body problems. Very few extensions to many-body systems 
have been considered so far. On the other hand, the primary 
motivation for studying low-dimensional systems is to under- 
stand the surprising propert ies of low dimensional many -body 
systems (see Refs. iITsI - UTIi for a review). Moreover, experi- 
ments are always performed in systems with a finite density; 
in many cases it may not be easy to disentangle many -body ef- 
fects from two-particle scattering. For example, confinement- 
induced molecules are relatively large on the BCS side of res- 
onance ifioll . Already for a modest density of fermions, dis- 
tances between particles may become comparable to the size 
of bound pairs, and the Pauli principle can have a strong effect 
on the collisional properties of atoms and, as a result, on the 
properties of CIR's. 

In this paper we provide a theoretical analysis of a many- 
body system composed of two-component fermions confined 
in two-dimensional geometries in the vicinity of a Feshbach 
resonance. We focus on quench-type experiments, where a 
noninteracting mixture is rapidly taken to the regime of strong 
interactions [1I8L I19i1 . We analyze many-body corrections to 



the energies of confinement-induced molecules and calculate 
the rate at which they are formed out of unbound atoms. 

One of the intriguing questions raised by recent experi- 
ments concerns the possibility of using fermionic systems 
close to Feshbach resonances for exploring many-body phe- 
nomena associated with strong repulsive interactions. For ex- 
ample, positive scattering length on the BEC side of the Fes- 
hbach resonance has been suggested as a route to observe the 
Stoner instability I.20i1 . While the first experiments by Jo et 
al. have been interpreted using a simple mean-field picture 
of such transition flBJlX], subsequent measurements showed 
that the system is strongly dominated by fast molecule forma- 
tion ll22ll . as predicted theoretically in Refs. lfT9ll23ll . In this 
paper, we demonstrate that dimensional confinement can have 
a dramatic effect on the dynamics of molecule formation. We 
find that the peak in the molecule formation rate should be 
shifted from the BEC to the BCS side of the resonance with 
increasing transverse confinement. Testing our predictions in 
experiments will help to di sting uish between different models 
of molecule formation |fl9l 1231] . 

Another conceptually intriguing aspect of the system we 
study is that one cannot use separation of energy scales to 
simplify the analysis. Typically when many-body systems of 
ultracold atoms are studied, it is assumed that one can start by 
solving a two-body problem to obtain the strength of contact 
interaction and then work with this contact interaction when 
analyzing the many-body problem. In our system the effec- 
tive two particle scattering can be strongly modified by the 
presence of other particles [19]. Hence, an accurate analysis 
of our system requires understanding of the interplay between 
few-body and many-body phenomena. 



II. VACUUM T-MATRIX 

Traditionally, two-body problems in low-dimensional ge- 
ometries have been analyzed using the Schrodinger equation, 
which can be simplified into two decoupled single-particle 
problems coiTesponding to the relative and the center-of-mass 
(cm.) motion |9, 10]. This approach is, however, difficult to 
generalize to the many-body case. In the presence of a filled 
Fermi sea, the cm. momentum of the scattering pair relative 
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to the Fermi sea is important and cannot be taken into account 
by a simple momentum boost. Therefore we re-examine the 
two-body problem in quasi-2D geometries by recasting the re- 
sults of Ref. 1 10] to the form of a T-matrix in vacuum. For a 
discussion regarding Feshbach resonances in low-dimensional 
systems, see Ref. ll24ll . We take the gas to be homogeneous in 
a 2D plane and assume a strong harmonic confinement in the 
transverse direction. 

We start from the full 3D scattering problem and use a 
contact interaction V5„j(r — r') = Vo(5(r — r') to describe 
the inter-particle interactions. In order to make the connec- 
tion to the many-body problem we do not separate relative 
and center-of-mass motion from the outset. This gives rise 
to a T-matrix which depends on energy hui as well as on 
the harmonic oscillator quantum numbers n ~ {^1,11,2) and 
n' = (^1,^2) corresponding to incoming and outgoing par- 
ticles. For the contact interaction the Lippmann-Schwinger 
equation takes a simple form 

Tlft'H = VnM' + V Vn,n"n^°]{u;) r«, ^,(0.), (1) 



where the polarization operator is given by 
dk 1 



" ^ ^ J {2n 



and we have denoted the trap frequency in transverse direction 
by u)^. The polarization operator in vacuum has the property 



n?''(LLi) — n^"''_|_„ (lu) fovn = (ri]^, 71-2 )■ We will utilize these 
two notations interchangeably when discussing the proper- 
ties of the many-body T-matrix. The dispersion is given by 
= h^k"^ /2m, and we measure energies and frequencies 
with respect to the zero-point energy hbJz- Thus, a confined 
particle in the lowest vibrational state with no in-plane mo- 
mentum is assumed to have zero energy. 

Since V-^^^ depends only on the relative motion of scattering 
particles, we write the matrix elements Vfi,n' in terms of the 
quantum numbers corresponding to relative (n^) and center- 
of-mass motion (iV) 



t(o) 



(3) 



Here Lp^ is a harmonic oscillator eigenfunction correspond- 
ing to relative motion and the harmonic oscillator length in 
the transverse direction is denoted by = ^Jh/moj^. The 
Clebsch-Gordan coefficients arising from the change of ba- 



sis are defined as C 
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Quantum num- 



bers ni,n2, Uj., and N are non-negative integers and energy 
conservation imposes condition Ui + n2 — N + n,. for the 
non-zero elements „ . 

The form of matrix elements Ki.n' suggests we look for a 
solution in the basis of relative and cm. quantum numbers 
and then go back to the original basis. We find that (see Ap- 
pendix A) 



N.i 



<^;jo)(^„, (0) 

X To{to- NloJ. 



(4) 



The structure of T2 ^, shows explicitly the decoupling of rel- 
ative and cm. motion. Furthermore, since the interaction po- 
tential depends only on the relative motion, the cm. quantum 
number does not change in the scattering and contributes only 
as shift to the energy of scattering particles. When the bare in- 
teraction Vq is eliminated. To is given by lfloll25ll (for details, 
see Appendix A) 



1 



where function w{z) is defined as 



[V2^iJa3jj + wiu}/uJz+iO+)], (5) 



w{z) — lim 



2W^ln4 



(2^-1)!! 
(2^)!! 



ln(£- z/2) 



(6) 

The double factorial is given by n!! = n • (n — 2) • (n — 4)..., 
and by definition (—1)!! = 0!! — 1. The two-body T-matrix 
has a series of poles corresponding to different values of the 
center-of-mass quantum number N. In particular, there is a 
bound state corresponding to iV = which exists for all a^jj 
and coincides with the Feshbach molecule deeply on the BEC 



side. Deeply on the BCS side of resonance 1 



'■3D 



< the 



energy of the confinement-induced two-body bound state has 
(2) a simple expression ~ —^hu^ q-V2tt ej\a.^^\ ^ where B = 
0.905 1 1]. In general the pole has to be computed numerically 
from Eq. (|5]l. 



III. MANY-BODY T-MATRIX AND COOPERON 

Let us next discuss the many-body effects in the forma- 
tion of confinement-induced molecules in quasi-2D geome- 
tries. The system is described by a many-body Hamiltonian 



k,n.(T 
1 

+ v 



(7) 



where ^ 



/iw^rtj and particles carry 2D 



momentum k as well as harmonic oscillator quantum number 
n^. We have also allowed a possible imbalance between the 
two fermion species. 

To incorporate the Pauli blocking to our analysis, we de- 
rive a T-matrix in the presence of Fermi sea (Cooperon). We 
approximate the full Bethe-Salpeter equation by taking into 
account the ladder diagrams and obtain 

Ti^|,(^, q) = Vn,n' + E Vn,n"'Rn-{^, q) T^^A^, q), 

n" 

(8) 

where we assume that the scattering particles can have finite 
cm. momentum q in the 2D plane. The full polarization op- 
erator n,,-j(a;, q) is of the form 
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dk 1 



^(4 



)-"f(4.«, l) 



(2^)2 fi,.(l + zO+)-4 



(9) 

Although cm. and relative motion become coupled in the 
presence of Fermi sea, we utilize insights from the two-body 
problem and look for a solution where cm. and relative mo- 
tion are at least partially decoupled. We find that the solution 
of Bethe-Salpeter equation can be written in terms of a T- 
matrix depending only on the cm. quantum numbers 



2^4 J2 qXc»',,„, <(o)^„,(o)r, 



N,N' 



N,N' 



(10) 

We use the two-body T-matrix Tq to renormalize the UV di- 
vergence associated with polarization operator (|9|l and obtain 
(for details, see Appendix B) 

(11) 

where the renormalized polarization operator is given by 



T^N,N' — ^ "ri,,,JV+ri,,-_ff C 



pn 

K,N+n^-K '^N,n 



X [ns(a.,q)-nf (^-^,)]. (12) 

The coefficients u„ ,„ are related to the zeros of the harmonic 
oscillator eigenfunctions [see Eq. ( IA4b 1 and they are given by 

1)!! 

II ^ 

n.m I — i r 

\Jn\ ml 

for even and non-negative n and m. Otherwise m„ , 

1 



(13) 



IS zero. 



q — 2^q '~PA ^PA' ^'^ order 



to correctly renormalize the UV divergence associated with 
the 2D momentum integral in Eq. (|9]l, we have to evaluate the 
two-body T-matrix such that the Fermi surface and finite cm. 
momentum are taken into account. This shifts the argument 
of Tq by in Eq. (fTTT i. 

Conservation of energy and parity impose selection rules 
for the allowed scattering processes and render matrix 1^ 
non-invertible. Since both T" and T^^^ share the same struc- 
ture as they also lack well-defined inverses and Eqs. (|4|i 
and ( [Tol l have to be solved in terms of matrices Tq and T 
which are both regular. The full solution retains all discrete 
energy levels in the transverse direction and although the most 
interesting 2D limit does not involve real processes via higher 
bands, virtual scattering processes become important near the 
Feshbach resonance. The general solution based on Eqs. (fTOli- 
(fTZt enables a systematic analysis of pairing instabilities from 
the strictly 2D regime at zero temperature to the confinement 
dominated 3D regime where temperature and Fermi energy 
become comparable with fiw^. 



IV. MOLECULE FORMATION 

To analyze the possible pairing instabilities, we assume that 
the system is initially spin balanced and weakly interacting. 



In the spirit of Ref. ifTsll . we consider an instantaneous quench 
where interactions are rapidly modified utilizing a 3D Fesh- 
bach resonance. The molecule formation is associated with 
the appearance of poles fiw fig + lA^ in the many-body 
T-matrix T^?,{Lu,q) [26]. We identify the real part fig as 
the binding energy of the molecule and the imaginary part 
as the growth rate of the instability toward formation of 
molecules 1 19]. 

Similarly to the 3D case ifioll . we find that the system ex- 
hibits an instability towards molecule formation via two-body 
processes as long as the Fermi sea can absorb the binding en- 
ergy of the molecules. This results in a sharp cutoff in the 
growth rate, see Fig. ^a). For a fixed t^/a^j^, the binding 
energy of molecules depends strongly on the ratio ep/huj^ 
and Fig. [Ttb) shows that the binding energy increases with 
increasing strength of the transverse confinement. The loca- 
tion of the peak value for the growth rate of instability can be 
varied by adjusting the ratio i^/a^j^ and, in particular, tight 
enough transverse confinement can move the pairing insta- 
bility completely to the BCS side. On the other hand, when 



Epjhw^ ~ 0.1 as in Refs. 012112711 . the pairing instability ex- 
tends to the BEC side and fast two-body processes dominate 
the three-body processes. When the molecule formation via 
two-body processes is no longer possible, the leading insta- 
bility is a three-body recombination which is suppressed for 
Fermi gases due to low densities and the Pauli principle. 
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Figure 1. (Color online) Growth rate of the pairing instability (a) 
and the binding energy of molecules (b) at zero temperature as a 
function of l^/a^o- The values of Sp/hu)^ are (from left to right) 
Ep/hjj^ = 0.0175, 0.025, 0.0375, 0.075, 0.2, and 0.4. 

The binding energy in vacuum is compared to the bind- 
ing energy at finite densities in Fig. |2] The relation be- 
tween vacuum and finite density binding energies depends 
again on £p/hjj^, and when the system becomes more three- 
dimensional (3D) (i.e. when Ep/hw^ increases), Pauli block- 
ing by the Fermi sea can result in a stronger binding of 
molecules. The crossover takes place roughly at ej;./Sw^ = 
0.5. In 3D gases many-body corrections always result in 
stronger binding [19.1 and Pauli blocking induced weaker 
binding is a manifestation of 2D physics. On the other hand, 
the binding energy of the molecules is larger than the binding 
energy of the Feshbach molecules existing on the BEC side of 
resonance. 

Finite temperature suppresses strongly the growth rate of 
pairing instability whereas the binding energy decreases more 
slowly with increasing temperature. In Fig. [3] the growth rate 
is shown at different temperatures for Ep/hw^ = 0.1 corre- 
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Figure 2. (Color online) Binding energy of molecules at zero tem- 
perature for Sp/I'ijjj^ — 0.025, 0.075, and 0.2 (from bottom to top). 
The binding energy in the presence of Fermi sea (solid lines) is al- 
ways smaller than the vacuum binding energy (dashed lines) for the 
parameters investigated. 

sponding to the experimental parameters of Refs. lfT2l l27ll . 
Pairing instability at the BCS side of the resonance is sensitive 
to the temperature since thermal fluctuations can easily break 
molecules at small binding energies. At high enough tem- 
peratures the pairing instability can become completely sup- 
pressed for weak attractive interactions. On the other hand, 
the cutoff in the growth rate at i^/a^o ~ 0.5 does not in 
general depend strongly on the temperature. We note that al- 
though the pairing instability can persist to quite high temper- 
atures, the critical temperature for the superfluid transition is 
typically much lower near the Feshbach resonance or deeply 
in the BEC regime [28,,29|]. 
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Figure 3. (Color online) The growth rate of the pairing instability as 
a function of i^/azD for ep/hw^ = 0.1. The temperature is (from 
top to bottom) T/Tp = 0, 0.1, 0.2, 0.3, 0.4, and 0.5. 

The vacuum T-matrix has several poles on the BEC side 
of the Feshbach resonance corresponding to the different cm. 
quantum numbers. This can in principle give rise to seve- 
ral pairing instabilities which show up as distinct poles in the 
Cooperon. For small Ep/fujj^ these poles do not coexist for 
given l^/a^jj since the Fermi sea is unable to absorb large 
binding energies. When Sp/huj^ increases, the poles start to 
overlap and additional poles with non-zero Ag become dis- 
cernible. However, these additional instabilities remain weak 



compared to the primary instability corresponding to the pole 
of the vacuum T-matrix with = 0. 

So far we have analyzed the pairing instability in the case 
of zero cm. momentum q. We find that the results remain 
qualitatively the same for finite q and in the spin-balanced 
case the most unstable mode is always at q = 0. However, 
the growth rate of instability decreases slowly as a function of 
|qr|; in a realistic quench experiment it is likely that molecules 
with a wide distribution of momenta are created. We find that 
finite q reduces the binding energy due to the smaller number 
of low-energy states that are available for scattering jsoll . In 
spin-imbalanced systems, the lowest energy state can shift to 
finite momentum llsTl [32ll . 



V. DISCUSSION 

We have studied pairing instabilities in spin balanced quasi- 
2D Fermi gases when interactions are dynamically quenched 
to the regime of strong interactions using 3D Feshbach reso- 
nances. We found that the pairing instability can be shifted 
to the BCS side of resonance by adjusting the axial confine- 
ment with respect to Fermi energy. Pauli blocking was found 
to renormalize significantly the binding energies, resulting in 
weaker binding in the 2D limit than that warranted by the two- 
body description. 

The growth rate of pairing instability can be measured by 
monitoring the atom loss |18] and the binding energy can 
be probed using rf spectroscopy |i3 HH [34ll . In a related 
work fssf], we argue that the recent experiment lfT2ll prob- 
ing the properties of 2D Fermi gases can be interpreted in 
terms of dynamically created polarons. Another recent ex- 
periment 1 34] measures directly the binding energies of the 
molecules and finds agreement with a theoretical prediction 
for the two-body bound states in ID optical lattices f36]. On 
the other hand, our calculation (Fig.|2]i predicts that the two- 
body bound-state energy should be significantly renormalized 
by the presence of the Fermi sea. The discrepancy could 
stem from the fact that our calculation probes an unpaired gas 
which is rapidly quenched to the strongly interacting regime, 
whereas in Ref. [34] the system corresponds to a strongly 
interacting gas in equilibrium with a large number of paired 
atoms. 

The T-matrix approach presented here can be used to 
probe the competition between polaron and molecule lIsTl 
[32i [37I [38I1 in quasi-2D systems and to investigate dimen- 
sional crossover from 2D to 3D |29, 34, 39, 40]. Our for- 
malism is also useful for studies of pair formation in other 
low-dimensional geometries and Bose gases. In particular, it 
can be used to investigate the effective three-body collisions 
induced by virtual excitations of the transverse modes Il4ll - 

El. 
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(b) follows from an analogous calculation. We note that the 
harmonic oscillator eigenfunctions in Eq. (IA2I) satisfy 



( (n-lVA f 

^ « I /^\i9 — rr—, for even n 
"^'^^ ^' 10, foroddn. 



Appendix A: Vacuum T-matrix 

We briefly discuss the technical details regarding the calcu- 
lation of the vacuum T-matrix as well as the many-body T- 
matrix. For simplicity, we set ?i = 1 in Appendices A and B. 

The form of the matrix elements in Eq. (O suggests the 
following ansatz for the vacuum T-matrix 



N.n ,n' 

(Al) 

Substituting Eq. (lAlb into Eq. ([T]l we obtain 



2^4 



E 



dk 



We observe that the cm. index comes only through the shift 
of energy. Thus we can take Tn{uj) —Tq{lo~ Nlu^) and for 
To{uj) we obtain 



1 V2^£^ 



E 



dk V2TT£jip„Mr 



(27r)2 ^ - - c^^n, + zO+ ■ 
(A2) 

We can calculate the integral in Eq. iA2i using the identity 

1 /•TOO 



A 



dTe^\ 



(A3) 



where =F = sgn(Re A). For Eq. (IA21 l we have two cases: (a) 
Re w < and (b) Re w > 0. We discuss case (a), and case 



Using the identity in Eq. ( IA3b . we obtain 

V- f dk V2^4|(^„^(0)p 



oo go 

dr > ^— — — ^ — e 
^ (271 !! 

n=0 ^ ' 



V47rr/ 



LO T \ 1/2 



Q \2sinna;,T/ V47rr/ 



The cm. part of the quasi-2D T-matrix satisfies, therefore, an 
equation 
1 



T^{lo) Vo Jo y 2sinhw^T V47rT/ ' 

(A5) 

The UV divergence associated with the original contact in- 
teraction is manifested as a singularity of the integrand in the 
limit r — >^ 0. We regularize this divergence using the 3D T- 
matrix, which is given by an analogous equation: 



1 1 r 

^3d(^) ^0 Jo V47rr/ 



(.+.0^) (n^y/^ ^^^^ 



We take w — > of Eq. ( IA6I ) to obtain 



m 1 



47ra3^ Vo 



dr 



(—) 



3/2 



(A7) 



Using this identity, we eliminate the bare interaction Vq from 
Eq. dASI l. This gives us an To{uj) which is manifestly free 
from UV divergences: 



Att I a^jj Jq X \ 2 sinh x J V 2a; / j 



The latter term in Eq. ( IA8I ) is the integral representation of the 
function w{lo/lu^ + iO+) defined in Eq. (|6]l. 



Appendix B: Many-body T-matrix 

To solve the Bethe-Salpeter equation (|8]l we generalize the 
ansatz in Eq. ( lAll i and assume that the many-body T-matrix 
is of the form 



I 

T^^l = V2^4 ^ C£^C^'.<<(0)^„;(0)T^,^,„ 

N,N' 

(Bl) 

where we have temporarily suppressed the frequency and mo- 
mentum arguments. The polarization operator satisfies the 
following useful identity: 
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n n' 



(B2) 



where — ^^k^Sp^^^F^i- Substituting the ansatz ( IB lb to 
the Bethe-Salpeter equation (O and using the above identity, 
we obtain an equation for the cm. part: 



N' 



^ K 



Vr 



^Y.^V)j,^j,Tk,n', (B3) 



2ttI 



K 



where matrices 2?o and T) are given by 

{'Do)n,K = 5n.K E "rv,n.<Vn,(^ " ^k). (B4) 



and Eq. (fTZt . respectively. Coefficients „ are given by 
Eq. ( fT3] l in the main text. Equation ( IB3I ) for T = (TJy ^/ ) can 
be written in a matrix form 



(B5) 



^ 1 V27r^^ _ _ 



Denoting T^°^ = diag(Tjv), where TNi^jj) — To{uj — Nuj), 
we obtain 



(B6) 



This gives us an equation for the many-body T-matrix such 
that the UV divergence associated with 2?o is renormalized 



^-1 ^ ^(o)-i _ ^ 



(B7) 



Equation ( |B7| ) is illustrated in more detail in Eqs. ( fTTI ) and ( fT2b 
of the main text. 
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